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A  method  is  described  for  computing  the  emissivity  1.) 

of  a  cylindrical ly  symmetric,  optically  thin  li^*ht  source, 
when  the  projected  intensity  profile  ■  ■JlF(f)  dy  is 

toiown  from  experiment.  The  unloiovm  function  is  taken  as 

a  series  expansion  in  terms  of  orthogonal  polynomials,  and  it 
is  shown  how  the  expansion  coefficients  can  be  determined 
U»h  This  procedure  yields  a  least-squares  smoothed  approxi¬ 
mation  f.Lr). 

Let  f(r)  be  the  emission  coefficient  (emission  per  unit  area)  of 


a  cylindrically  symmetric,  optically  thin  light  source  of  unit  radius. 
Side-on  measurements  of  the  emission  give  the  projected  intensity  pro^ 
file  l(x),  which  is  related  to  f(r)  oy 


I(x) 


r 

■J/- 


f(r)  dy  «  2 


1 


The  inversion  of  oq.  (l),  giving  f(r)  in  ter*ns  of  l(x),  is 


I 


l(y)  X  ex 
(x^-r^)^ 


The  i  roblc::.  of  dot .^rrin' 'p:;  f(r)  "(x;  kr.ovm  frox.  ex;  riru<  n 

hio  been  routed  by  .-ov'cral  'authors,  ^  ^  mostly  by  seme  method  of  stri: 

.  2 

integration,  .^n  '••.•cur-ite  nurr...r,cil  m  )tho  I  h  -.n  .‘‘;ivon  by  Bockast-  n, 


/■ 


2 


The  x-axis  is  divided  in  a  n\mbor  of  equal  int. rivals,  :ind  a  third-doGrec 
polynomial  is  fittsd  to  the  l(x)  curve  within  each  interval.  E*..  (2) 
is  used,  and  a  matrix  is  obtained  that  transforms  a  set  of  values 
of  I  to  values  of  f, 

f(r^)  “21 

Bockasten  also  discusses  the  influence  of  rc^ndonl  ax^rors  in  tho  values 
of  I(x). 

The  method  described  in  this  note  should  be  advantageous  wheneve}’ 
the  l(x)  curve  has  a  more  or  less  irregular  shape,  and  a  properly  smoothed 
approximation  to  f(r)  is  desired.  The  urJaiovm  function  f(r)  is  e::par.dod 
in  a  series  of  orthogonal  polynomials.  The  expansion  coofficientn  can  ho 
determined  from  l(x),  and  in  this  way  a  least-squares  sj’nocthed  a:  ro'O- 
mation  to  f(r)  can  bo  found  directly. 

A  few  general  consequences  of  eqs.  (l)  :and  (2)  arj  of  interest. 
Obviously  l(x)  is  an  even  furiction,  and  consider  the  interval 
0^1.  only.  The  behavior  of  f(r)  near  r«=^  1  io  very  sensitive  lo  ti  e 
behavior  of  l(x)  near  x  ^  1  .  If  l(x)  Is  proi  ^rtio’^al  to  near 

x»1,  then  f(r)  will  bs  proportional  to  (l-r)°*"  near  r^*1.  IV: 
instance,  b  *  0./9.  1  ^  0.5>  ^^^d  b  =*  0,5‘  cerro spend  to  iiifintts  ,  r  .. 
and  zero  valu.^s  of  f ( 1 )  s  rospec tivcly .  Further,  it  is  s.'on  fre*..  (2) 
that  1  discontinuity  r.n  l/(x)  at  x-^0  implies  an  Infini'.te  value  of  ^^O). 

V/o  will  assume  that  functions  l(x}  'jr.a  f(r)  arc  finite,  -o  that 

I"(C)  ..  0,  -nd  I;1)  «  0  (wi^h 

Fc  now  choose  to  fit  in  r.  polyi-o'-ial  ieyree  h  ,  tv 

in  such  a  way  that  the  error  squared,  into^Ti’ated  over  the  croa'^-  v.notio.i 
r^  1 ,  ir.  -’’'hi::  tins. 


3 


I 


r  dr  ■  min.,  M  fixed. 


The  orthogonal  polynomials  suitable  for  this  purpose  are  tha  Lof^endro 


A  P 

polynomials  ^  P^(t),  T-rith  t  ■  2r  -  1  .  The  first  few  of  those  axe 
-  1,  P^  -  2r^-1  ,  Pg  -  6r^-6r^  +  1,  P^  -  20r^-  30r^  +  12r^-1. 


They  satisfy  the  orthogonality  relation 


2( 2m+l ) 


(2r2-  1)  P„(2r2-  1  )  r  dr  .  ^  , 


BO  that  the  series  expansion  of  f(r)  is 


f(r)  -  H  p 


1 

•  2(2m+l)  ^  f(r)  P^(2r^-1)  r  dr. 


If  the  series  (6)  is  terminated  at  in«M,  the  result  is  a  2i4-do^oc 
polyn<»nial  f^,  which  satisfies  the  condition  (4)# 

We  substitute  the  series  (6)  in  oq.  (l),  integrate  tom  by  terti., 


and  obtain 


l(x)  «  2(2in+l)  a^sin  r  ( 2in+1  )c03  ^x  |, 


or,  putting  x  •  cos9, 


l( cos©) 


2(  2in+1 )  '  2m+1 


'  The  functions  U  (x)  ?  sin(n  cos  ^x)  are  related  to  tho  Choby3i.-j  . 

^  “"I 

polynomials  T  (x)  a  C03(n  cos  x).  They  satisfy  the  orthogonality  r  ' 
^  1  ^2  V 

tion  \  U,  U  (l-x*^)  ~  dx  ■  0,  Ic/n,  and  vanish  at  the  points  x  =»  +  ', 
X  n'  '  4 

U  (x)  is  an  even  function  if  n  is  odd.^ 
n 


4. 


It  can  be  shown  from  eqs.  (l)  and  (t)  that 

p2 

2(201+1)"''  a  «  (4/F)  )  I(oo30)  3in(2m+l)9  d«,  (10) 

80  that  eq.  (l)  (with  x  "■  cosO)  transforms  the  Le/jendro  aariar  of  f(:.*} 
into  the  Fourier  series  of  I(co3Q)  (odd  sine  terms  only).  It  may  be  notoo. 
that  the  total  omission  from  the  light  source  is 

21f  (  f(r)  r  dr  =  2  ^  l(x)  dx 

0  0 

It  i3  also  seen  that,  if  l‘jj(r)  satisfies  the  condition  (4)»  then  the 
corresponding  function  (the  series  (8)  terminated  at  m»I>i)  will 

be  the  "best’ approximation  for  l(x),  in  the  sense  that  the  integral 
1 

J  (I  (1  -  dx 

0 

is  ninimizod. 

5 

The  standard  methods  of  Fourier  analysis  can  nov;  be  used  fo^ 
determining  the  coefficients  One  possible  pr' cedure  is  is  follow:: 

The  series  (9)  is  torminatod  at  some  m  i 

...1 

Ijj(coE^)  e*  2(2n+l)  r. '.n(  ?;n+1  )G,  (1?) 

ivr  0 

and  is  recuiroi  to  ccincido  ;rlth  I  in  thcs>>  points  in  interval 

i. 

0-^9^. ^/2  vrhore  the  ;'irst  noglegto to:m  ( :jir/ 211+3)9)  vanishe.i.  Those 

points  are 

'  -  ;'J^/(2J+3}.  'X  -I  ,  2-.-.  .'+1. 


ind  the  0?.  ali  Ie 


5. 


a  •*2  ;S  l(oos  ■  )  sin  1...  K  .  (13) 

“  aJ+3  k-1  2N+3  2N+3 

This  result  can  also  he  obtained  by  evaluating  the  integral  (10)  by  the 
trapezoidal  rule.  Sq.  (13)  i*  exact  if  f(r)  Is  an  even  polynomial  of 
degree  2(N+1 )  or  less,  so  that  the  series  (6),  (8),  and  (9)  contain  no 
terms  beyond  m  ■  N+1 , 

With  the  a  given  by  oq.  (13),  the  2M-degree  polynomial  approxima- 
n 

tion  for  f(r)  becomes 

V-)  -  H  \  Pj2r2-1)  .  ^  Vr)  Kcos  ^  )  (I4) 

m-0  k-1 

vrhere 

"  ailT  i  sin  (15) 

m-o 


As  an  example  we  consider  the  fionction 


I  »  ( 1  -  x^)^  -  sin'^©, 

(16) 

corresponding  to 

f(r)  -  (8/3W)(l-r2)^/2 

(17) 

The  coefficients  a  f  computed  from  oq.  (l3)j  ^ro  listed  in  Table  !•  The 

Hi 

oxict  a^,  given  by  eq.  (j)  or  (lO)  are  entered  in  the  last  line.  They  ar^ 

^m  "  'V  ( 2r.-i)(  2m--! )( 2in  <-3  }(  2m+57  ’ 

showing  fairly  rapid  convergence.  The  convergence  slowest  neir  the 
points  and  r-1 ,  where  jp^j  1 ,  A  seven-point  -analysis  (n--6)  is 

aufficient  to  reduce  the  error  in  f(r)  to  about  0.0006  for  r«0  and  0.001 
for  r« 1 . 


Table  2  shows  the  same  for 


6 


I  ■  1-x^  ■  sin ^6, 

(19) 

oorresponding  to 

f(r)  -  (2/fr)  (1-r^)^. 

(20) 

The  exact  coefficients  a  are 

ID 

^  4  , _ 1 _ 

(21) 

*m  W  (2m-1)(2m+3)  ’ 

so  that  the  convergence  is  fairlj"  slow.  Clearly,  the  behavior  of  the 
function  (20)  near  r»1  cannot  bo  rerresontod  accuratoly  by  a  low-dogroo 
polynomial.  However,  a  seven-point  suialysia  (N*6)  suffices  to  give  an 
error  less  than  about  0.002  for  rii,0,96, 

Bxperimental  measurements  will  commonly  yield  l(x)  curves  with 
irregular  fluctuations,  rather  than  smooth  functions  of  the  type  (16) 
or  (l9)»  In  such  cases  it  should  bo  an  advantage  of  this  method  that, 
with  a  suitable  choice  of  N  and  M,  a  properly  smoothed  f(r)  is  obtained 
directly. 

The  influence  of  random  errors  in  the  values  of  l(x)  can  be  found 
as  follows.  If  the  error  (standard  deviation)  of  I  at  each  point  is/}I, 
eq.  (13)  gives  for  the  error  Aa^^^  of  a^, 


Biti^  (ai) 


2 

f 


or 


^  2n+1 

- - Al. 

(2Ii+3)^ 


(22) 


For  the  fimction  (19),  asstuning  «  0.002  and  11  *  6,  wo  see  from  cq.  (22^ 

and  Table  2  that  the  error  in  a-  is  comparablo  to  itself,  so  that  this 

5  5 

soiiroe  of  error  is  more  important  than  the  fairly  slow  convergence  of  the 
Legsndro  series. 

Prom  ©q.  (14),  the  error  in  f,.  is  ^;ivon  by 


7 


which  cannot  ho  oxpressed  in  a  simplo  form.  Wo  can  obtain  an  ostimato 
for  Af||  from  tho  expression 

m«o 

which  would  be  correct  if  thsAa^j^  woro  uncorrolated. Afjj  will  bo 
largest  at  the  points  r«0  eind  r»1 ,  whore  jP^|  ■  1 .  We  find  from 
eqs.  (22)  and  (24)> 

For  intorrnodiate  valuos  of  tho  error  will  bo  smaller,  although  no 

2 

simple  formula  can  be  givens  For  instance,  for  r  «  0s5  K  «  6, 
eqs.  (22)  and  (24)  yield /if g  -  5.9  x  (2N+3)”^I,  while  tho  end-point 
value  given  by  oq.  (25)  is  21  X  (2N+3)"W  .  For  M  -  10  the  factor  is 
9.0  and  42,  respsotively. 

Aoknowledffcont. 

The  author  is  indebted  to  Dr.  K.  Bockaston  for  valuable  discussions 
and  comment 8. 
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Table  1 

Expansion  ooofficionts  of  the  function  (8/3f)(l  -  ,  approximate 

values  from  eq.  (13)  with  N  -  0,  1,  3,  6,  and  (last  line)  correct  values 
from  oq.  ( 18), 


N 

^0 

^2 

^3 

^4 

"5 

^6 

0 

0.32476 

1 

0.33931 

-0.44089 

3 

0.33957 

-0.43660 

0.08062 

0.00951 

6 

0.33953 

-0.43657 

0.08083 

0.01027 

0.00306 

0.00117 

0.00046 

(o^ 

0.33953 

-0.43654 

0.08084 

0.01029 

0.00305 

0.00124 

0.00061 

Table  2 

Expansion  coefficionts  of  the  function  (2/lf)(l  -  r  approximate 

values  from  eq,  (13)  vdth  N  «  0,  1,  3»  6,  and  (last  lino)  correct  values 
from  eq.  ( 21 ) . 


V 

^0 

^2 

^^3 

^4 

"5 

^6 

0 

0.4330 

1 

0.4253 

-0.2437 

3 

0.4245 

-0.2538 

-0.0580 

-0^0216 

6 

0.4244 

-0.2546 

-C.O6O3 

-0.0277 

-0.0154 

-0.0089 

-C.OO44 

(0^ 

0.4244 

-0.2546 

-r.0606 

-0.02? 3 

-C.OI65 

-0.0109 

-0.0077 

